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Abstract: Nucleotide sequences of segments of the mitochondrial control regions were analyzed to infer the phyloge- 
netic relationships among 7 macaques. High nucleotide diversity in Macaca assamensis and relatively low diversity in 
M . thibetana were found. Based on the ML tree from control regions, species in our study can roughly be sorted into 
three species groups except for the phylogenetic position of M . fascicularis , i.e., silenus group, including M . leonina ; 
sinica group, including M .arctoides, M.assamensis, and M.thibetana; and fascicularis group, including 
M . mulatta and M .cyclopis. A discrepancy between earlier studies (Fooden & Lanyon, 1989; Tosi et al, 2003a; 
Deinard & Smith, 2001; Evans et al, 1999; Hayasaka et al, 1996; Morales & Melnick, 1998), our result supported 
the hypothesis that M . fascicularis diverged earlier than M . leonina . Mitochondrial paraphyly in eastern M . mulatta (with 
respect to M . cyclopis) and eastern M . assamensis (with respect to M . thibetana ) were clearly observed in our study. In 
accordance with the results of Y chromosome, allozyme, nuclear genes and some morphological data (Delson, 1980; 
Fooden & Lanyon, 1989; Fooden, 1990; Tosi et al, 2000, 2003a, b; Deinard & Smith, 2001), our study on con- 
trol region sequences supported M . arctoides to be classified into the sinica group. However, this result disagreed with 
the previous mtDNA studies (Hayasaka et al, 1996; Morales & Melnick, 1998; Tosi et al, 2003a). 
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基于 线粒体 控制 区 序列 的 猕猴 属 系统 发 育 研究 
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2. 云南 大 学 生物 资源 保护 与 利用 重点 实验 室 ， 云 南 昆明 650091; 3. 中 国 科学 院 研究 生 院 ， 北 京 100039) 


摘要 : 通过 线粒体 部 分 控制 区 DNA 序列 数据 探讨 7 种 猕猴 属 物种 的 分 子 系统 发 育 关系 。 结 果 表 明 熊 猴 的 
核 苷 酸 多 样 度 最 高 ， 而 藏 鹊 猴 核 苷 酸 多 样 度 较 低 。 基 于 控制 区 序列 数据 所 构建 的 最 大 似 然 树 ， 不 考虑 食 德 猴 的 
位 置 ，7 种 猕猴 物种 可 粗略 地 分 为 3 个 种 组 ， 即 狮 尾 猴 组 〈 包 括 北平 顶 猴 ) 、 头 中 猴 组 〈 包 括 红 面 钦 、 熊 钦 和 藏 
PER) MARERA (AAA RA SBR). SHA (Fooden & Lanyon, 1989; Tosi et al, 2003a; Deinard & 
Smith, 2001; Evans et al, 1999; Hayasaka et al, 1996; Morales & Melnick, 1998) 的 结果 不 同 ， 我 们 的 结果 支 
RRR CIP RROCEHER, RARER (相对 于 台湾 猴 ) 和 东部 熊 猴 FRR) 出 现 并 系 。 与 
Y 染色 体 、 等 位 酶 、 核 基因 以 及 部 分 形态 学 数据 推测 的 结果 (Delson, 1980; Fooden & Lanyon, 1989; Fooden, 
1990; Tosi et al, 2000, 2003a, b; Deinard & Smith, 2001) 一 致 ， 红 面 猴 应 归于 头巾 猴 组 ， 但 此 结论 与 前 人 
(Hayasaka et al, 1996; Morales & Melnick, 1998; Tosi et al, 2003a) 依据 线粒体 得 到 的 结果 有 较 大 分 歧 。 
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The genus Macaca is the most widely distributed 
non-human primates with its 19 extent species extend- 
ing from northwestern Africa to a number of islands and 
continental areas in southern and eastern Asia 
(Groves, 2001; Abegg & Thierry, 2002). 

On the basis of male genitalia, Fooden (1976) 
classified the macaques into four species groups: 
silenus- sylvanus group, including M . sylvanus , M . sil- 
enus, M .nemestrina and the Sulawesi macaques; 
group, 
M . assamensis, and M .thibetana; arctoides group, 
including M .arctoides ; and fascicularis group, in- 
cluding M .fuscata, M .mulatta, M .cyclopis, and 
M . fascicularis. However, Delson (1980) established 
his classification system with a little difference: 
M .arctoides falled within the sinica group and 
M . sylvanus formed its own group. The viewpoint had 


sinica including M.sinica,- M .radiata , 


been supported by many followers (Nozawa et al, 
1977; Cronin et al, 1980; Melnick & Kidd, 1985; 
Fooden & Lanyon, 1989; Tosi et al, 2000, 2003a; 
Deinard & Smith, 2001). However, Groves (2001) 
divided the genus into six species-groups based on mor- 
phologic traits: 1) M.sylvanus group: monotypic; 
2) M . nemestrina group: M . nemestrina , M . leo- 
nina, M.silenus and M.pagensis; 3) Sulawesi 
group; 4) M .fascicularis group: M .fascicularis and 
M .arctoides; 5) M.mulatta group: M .fuscata , 
M . mulatta and M . cyclopis and 6) M . sinica group: 
M . sinica , M . radiata , M . assamensis , and M . thi- 
betana . Zhang & Shi (1993) recovered Fooden's clas- 
sification with the help of mtDNA restriction analyses. 
Furthermore, the result of Morales & Melnick (1998) 
was less concordant with the separation of the species 
into four distinct species groups. In the tree based on 
mtDNA sequences, the bootstrap values of some clades 
were too low ( < 5096) (Hayasaka et al, 1996; Tosi 
et al, 2003a). Therefore, further investigation and 
more data of mtDNA will be useful to resolve these 
problems. In the present study, we analyzed segments 
of the control region sequences to obtain further in- 
sights into the phylogenic relationships of genus Maca- 
ca. 


1 Materials and Methods 


1.1 Materials 
thirty-eight individuals representing 7 species: 
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M .cyclopis, M.leonina, M . thibetana, M .assame- 
nsis, M .arctoides, M.mulatta, and M. fascicularis 
(Table 1) were investigated. Baboon was served as 
the outgroup. Total genomic DNA was extracted from 
whole blood or other tissues using the standard phenol/ 
chloroform protocols. 

1.2 MtDNA sequencing 

The mtDNA control region was amplified using the 
PCR primer LqqF ( TCCTAGGGCAATCAGAAAGAA - 
AG) and TDKD ( CCTGAAGTAGGAACCAGATG ) . 
Each mixture contained at least 25 ng DNA template, 
5 uL 10 x reaction buffer, 1 pmol/L forward and re- 
verse primer, 2 mmol/L mixed DNTP, 3 mmol/L Mg- 
Ch, 0.25 unit of Taq polymerase enzyme, adding 
sterile distilled water to 50 uL and overlaid with miner- 
aloil. Amplification was performed on a RoboCycler 
Gradient 40 (Stratagene) thermal cycler. 

The cycling conditions for all markers began with 
an initial 2 min denaturation at 95 C; followed by 35 
cycles of 1 min denaturation at 94 C, 1 min annealing 
at 60 % and 1 min of extension at 72 C; a 5 min fi- 
nal extension at 72 % concluded each reaction. 

The PCR products were tested using agarose gel 
electrophoresis and purified with the Gel Extraction Mi- 
ni Kit of Warson Co. (Shanghai, China) following the 
manufacturer's instructions. Purified PCR products 
were labeled using the BigDyeTM Terminator cycle se- 
quencing kits (PE Biosystems, Foster City, CA) with 
an ABI 377DNA automatic sequencer using the same 
primers as used in the PCR reactions. Cycle sequenc- 
ing reactions were performed following the instructions 
provided by the manufacturer. Both forward and re- 
verse directions were sequenced. 

DNA sequences were edited using DNASTAR 
(DNASTAR Inc.) and aligned by CLUSTAL W pro- 
gram (Thompson et al, 1994). 

1.3 Data analyses 

Base compositional information and the nucleotide 
diversity (x value) for control region were estimated by 
MEGA 2.1 (Kumar et al, 2001). TREE-PUZZLE 
(Strimmer & von Haeseler, 1996) was used to con- 
struct a. maximum-likelihood (ML) tree for the se- 
quence data. An optimal model that fit the data for ML 
trees was determined by using Modeltest 3.0 ( Posada 
& Crandall, 1998). HKY + G was the most appro- 
priate for ML analysis . Settings for the HKY + G 
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Table 1 Samples used in this study and nucleotide diversity within each taxon 









































No. Samples Origin Origin abbr. GenBank no. Samples size x value! (96) 
1 M . cyclopis Taiwan TW AY682594 1 关 
2 M. arctoides Yunnan, China YN AY682588 6 4.1 
3 M . arctoides Yunnan, China YN AY 682589 
4 M . arctoides Vietnam V AY682591 
5 M . arctoides Vietnam V AY682592 
6 M . arctoides South of China SO AY682590 
7 M . arctoides Yunnan, China YN AY682593 
8 M . assamensis Vietnam V AY682611 12 5.9 
9 M . assamensis Vietnam V AY682613 
10 M . assamensis Vietnam V AY682617 
11 M . assamensis Yunnan, China YN AY682614 
12 M . assamensis Yunnan, China YN AY682615 
13 M . assamensis Yunnan, China YN AY682616 
14 M . assamensis Yunnan, China YN AY682612 
15 M . assamensis Myanmar MY AY682619 
16 M . assamensis Myanmar MY AY682621 
17 M . assamensis Myanmar MY AY682620 
18 M . assamensis Myanmar MY AY682622 
19 M . assamensis South of China SO AY682618 
20 M . fascicularis Vietnam V AY682595 4 3.4 
21 M . fascicularis Vietnam V AY682596 
22 M . fascicularis Vietnam V AY682597 
23 M . fascicularis Vietnam V AY682598 

24 M . thibetana Sichuan, China SC AY682610 5 0.6 
25 M . thibetana Sichuan, China SC AY682607 

26 M . thibetana Sichuan, China SC AY682608 

27 M . thibetana Sichuan, China SC AY682609 

28 M . thibetana Sichuan, China SC AY682606 

29 M . leonina South of China SO AY682623 4 5.6 
30 M . leonina South of China SO AY682624 

31 M . leonina Yunnan, China YN AY682625 

32 M . leonina Myanmar NY AY682626 

33 M . mulatta Fujian, China FJ AY682600 6 4.5 
34 M . mulatta Fujian, China FJ AY682599 

35 M . mulatta Fujian, China FJ AY682601 

36 M . mulatta Fujian, China FJ AY682602 

37 M . mulatta Yunnan, China YN AY682603 

38 M . mulatta Guizhou, China GZ AY682604 

39 Papio Africa AF AY682605 1 关 


1 Kimura-2-parameter model was used for mtDNA data analysis; * Species with one sample was omitted in the present study. 


model were as follows: 


- InL = 3428.3647; base fre- 


quencies = (A 20.3253, C 2 0.2825, G = 0.1222, 
11.9858, and the 
shape parameter of the gamma distribution = 0.2221. 


and T = 0.2699), Ti/tv ratio 


2 Results and Discussion 


Thirty-nine individuals were successfully se- 
quenced and 33 haplotypes were identified. Sequences 
had been deposited into GenBank (Table 1). The 
strong bias against guanine was observed in our ampli- 
fied sequences (G = 12.8%, A = 31.1%, C = 
28.7%, T = 27.4%), which characterized mtDNA 


rather than Numts (the nuclear insertions of mtDNA se- 
quences). In the 541 bp sequences of the control re- 
gions, a total of 230 variable sites and 174 parsimony 
informative sites had been found. The nucleotide diver- 
sity (x value) of control region in M . assamensis was 
slightly higher than in any other species, however, 
M . thibetana showed a relatively low diversity (Table 
1). 

As illustrated by the ML tree based on control re- 
gions ( Fig.l ), except for M.assamensis and 
M .mulatta, the different haplotypes found within a 
single species generally clustered together on a branch 
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65 M.mulatta (37;YN) 
0.1 M. mulatta (38;GZ) 


Fig.1 The maximum-likelihood (ML) tree based on control region sequences used by TREE-PUZZLE 
Only reliability percentages (after 1 000 quarter puzzling steps) above 50% are indicated. Numbers and letters in parentheses corre- 
spond to sample identifications and origin abbreviation outlined in Table 1. 


of the tree that did not include any haplotypes from oth- 
er species. Contrasting with the classification system of 
Delson (1980), species examined could roughly be 
divided into three species group except for the phyloge- 
netic position of M . fascicularis: silenus group, in- 
cluding M . leonina ; sinica group , including M . arc- 
toides, M .assamensis, and M .thibetana; and fas- 
cicularis group, M . mulatta and 
M .cyclopis. A discrepancy with earlier studies (Food- 
en & Lanyon, 1989; Tosi et al, 2003a; Deinard & 
Smith, 2001; Evans et al, 1999; Hayasaka et al, 
1996; Morales & Melnick, 1998), M .fascicularis 
diverged early in our ML tree based on control region 
with low reliability support (70%, Fig.1). 

Two macaque species were paraphyletic in this 


including 


study: M. mulatta (with respect to M .cyclopis) and 
M . assamensis (with respect to M . thibetana ). Previ- 
ous mtDNA studies had suggested a paraphyletic status 
for particular macaque 


species, for example, 


M .assamensis ( Hoelzer et al, 1992), M. mulatta 
(Melnick et al, 1993; Hayasaka et al, 1996; 
Morales & Melnick, 1998; Tosi et al, 2002, 
2003a) and M .nemestrina (Evans et al, 1999). In 
the paraphyly of M.mulatta (with respect to 
M.cyclopis) and M.assamensis ( with respect to 
M . thibetana ), earlier studies ( Hoelzer et al, 1992; 
Melnick et al, 1993; Hayasaka et al, 1996; Morales 
& Melnick, 1998; Tosi et al, 2002, 2003b) sug- 
gested a major separation between the eastern and the 
western lineages of M . mulatta and M . assamensis oc- 
curring along the Bramaputra river valley. The eastern 
M.mulatta lineage was more closely related to the 
M.cyclopis than to the 
M .assamensis was composed of eastern and western 


western lineage. 
sub-species: M.a.assamensis and M .a .pelops, re- 
spectively ( Fooden, 1982). Hoelzer et al (1992) 
documented that the eastern M.assamensis lineage 
(M.a.assamensis ) was more closely related to 
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M thibetana than to the western lineage (reviewed in 
Hoelzer & Melnick, 1996). 

Our results, however, was entirely out of expec- 
tation. In this study, we just had the eastern samples 


of M.mulatta and M.assamensis. M. mulatta and 


M .assamensis were still paraphyletic with respect to 
M .cyclopis and M . thibetana. M . mulatta cluster 
(including M . cyclopis ) clearly diverged into two main 
clades. M.mulatta of Fujian origin (eastern) di- 
verged first and formed one clade. The other clade 
contained M . cyclopis and M . mulatta of Yunnan and 
Guizhou origins (eastern). The departure from intra- 
specific monophyly was seen at its most extreme in the 
case of M . assamensis. Twelve individuals were suc- 
cessfully sequenced and three main clades have been 
occurred among them. The samples collected from the 
same place did not define an exclusive clade. Three 
M . thibetana haplotypes derived from five samples 
grouped an independent branch and swamped into the 
M . assamensis haplotypes. We could not continue fur- 
ther addressing about the paraphyly mechanism of 
M . mulatta and M . assamensis because we were short 
of the western M.mulatta, western M . assamensis 
(M.a.pelops) lineages and other species within the 
fascicularis group and the sinica group. Additional 
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本 书 论述 了 分 子 生 物 学 与 植物 分 子 生 物 学 简 史 、 植 物 的 基因 组 与 基因 组 学 、 基 因 的 结构 、 基 因 的 表达 





与 调控 、 基 因 的 功能 分 析 、 植 物 细胞 的 信号 转 导 、 植 物 细胞 的 分 化 和 器 官 发 育 、 植 物 生物 技术 等 内 容 ， 并 
附 有 植物 学 和 分 子 生 物 学 词汇 。 所 有 各 章 均 先 介绍 基本 知识 和 理论 ， 然 后 讲述 相关 的 专题 。 这 些 专 题 都 是 
近年 来 植物 分 子 生 物 学 的 重大 突破 ， 因 此 具有 很 高 的 科学 意义 和 现实 意义 。 同 时 ， 通 过 讲评 经 典 的 科学 论 
文 来 达到 能 力 培养 的 目的 。 本 书 图 文 并 茂 ， 适 合作 为 大 学 高 年 级 或 研究 生 的 教学 参考 书 。 
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数量 生态 学 的 内 容 十 分 丰富 。 本 书 内 容 涉及 到 数量 生态 学 的 各 个 方面 ， 包 括 过 去 半 个 世纪 来 有 关 数 量 
生态 学 的 主要 发 展 ， 全 面 介 绍 现代 数量 生态 学 理论 和 方法 ， 特 别 注 重 介绍 国际 上 最 常用 的 和 近年 来 发 展 的 
新 方法 。 人 全书 反映 了 数量 生态 学 发 展 的 概貌 ， 同 时 也 展现 了 我 国 在 这 一 领域 的 研究 成 果 。 对 于 有 关 学 科 的 
科技 工作 者 和 研究 生 ， 本 书 的 内 容 已 基本 能 够 满足 需要 。 书 中 所 讲 的 方法 ， 除 植物 生态 学 外 ， 可 以 普遍 应 
用 于 生物 学 、 农 学 、 林 学 、 地 理学 、 医 药 、 环 境 科 学 等 学 科 领 域 。 
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